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Abstract. We renormalize to three loops a version of the Thirring model where the fermion 
^ ■ fields not only lie in the fundamental representation of a non-abelian colour group SU(N C ) but 

also depend on the number of flavours, Nf. The model is not multiplicatively renormalizable 
in dimensional regularization due to the generation of evanescent operators which emerge at 
each loop order. Their effect in the construction of the true wave function, mass and coupling 
constant renormalization constants is handled by considering the projection technique to a new 
order. Having constructed the MS renormalization group functions we consider other massless 
independent renormalization schemes to ensure that the renormalization is consistent with the 
equivalence of the non-abelian Thirring model with other models with a four-fermi interaction. 
One feature to emerge from the computation is the establishment of the fact that the SU(Nf) 
Gross Neveu model is not multiplicatively renormalizable in dimensional regularization. An 
evanescent operator arises first at three loops and we determine its associated renormalization 
constant explicitly. 
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1 Introduction. 



Models with a four-fermi interaction have played an important role in exploring fundamental 
properties of physical quantum field theories. For example, the two dimensional O(N) Gross 
Neveu model, is similar to QCD in that it is asymptotically free. Another two dimensional 
four-fermi theory which shares similar properties to the Gross Neveu model is the Thirring model 
with a non-abelian symmetry, Q. It too is asymptotically free and additionally is more closely 
related to QCD itself. For instance, when the fermion fields of the non-abelian Thirring model, 
(NATM), have the same quantum numbers as the quarks of QCD, with Nf flavour dependence 
in addition to being in the fundamental representation of the colour group, it has been argued 
that it is equivalent to QCD in the large Nf limit, ||. More precisely, Hasenfratz and Hasenfratz 
demonstrated that by taking the NATM in its d-dimensional extension then integrating over the 
leading quark loops which dominate when Nf is large, the usual three and four gluon vertices 
of QCD are recovered as one approaches four dimensions. Moreover, graphs with a quark loop 
and five or more external gluon legs were shown to be irrelevant in the four dimensional limit 
when Nf was large, || . This equivalence can be understood better in the context of the critical 
point renormalization group. In effect one has a fixed point equivalence between QCD and the 
NATM at the d-dimensional Wilson-Fisher fixed point of QCD which is completely analogous 
to the d-dimensional equivalence at the Wilson-Fisher fixed point of O(N) </> 4 theory and the 
two dimensional 0(N) a model both treated in ^-dimensions. This underlying <i-dimensional 
property of QCD has been exploited in the large Nf expansion to compute 0(1/Nf) corrections to 
various renormalization group functions in the NATM in <i-dimensions, || |(|. More recently, 
results have also been determined at 0(1/Nf), J7j. As the NATM is equivalent to QCD then 
the information contained in these scheme independent critical exponents also relate to the 
perturbation theory of QCD. The calculational advantage of using the NATM is the absence 
of the graphs with gluon self interactions which substantially reduces the number of Feynman 
diagrams which need to be computed. Given this intriguing connection between the NATM and 
QCD it is important to understand the quantum and renormalization properties of the former 
theory in more detail. For instance, in QCD the four loop /3-function and mass anomalous 
dimensions have recently been determined in MS in |§, These built on the earlier loop 
calculations of [JO], 11, 13]. However, the situation in the NATM is unfortunately not as 



advanced. The model has only been renormalized to two loops, [14, 15]. Therefore, we have 
undertaken to renormalize the NATM at three loops which is the main topic of this article. 
Moreover, it is worth emphasising that aside from the QCD connection the two dimensional 
NATM deserves study in its own right given its rich structure and connection with integrability. 
Indeed the computation of 5-matrix elements is closely connected with the latter property and 
their computation in perturbation theory can only proceed when the model has been rendered 
finite. 

As will become clear in the discussion, however, performing the renormalization is a highly 
non-trivial exercise. Unlike the majority of useful quantum field theories the NATM is not 
multiplicatively renormalizable in dimensional regularization, ]jl5|| . Though the original two 
loop calculation of [|14|] was performed in strictly two dimensions using a cutoff regularization, 
it does not seem appropriate to us for technical reasons to develop that regularization for a 
three loop calculation. Further, dimensional regularization is widely used and has been applied 
to the renormalization of four-fermi theories before at three loops in (l^, [l?], 19]. The 



lack of multiplicative renormalizability is manifested by the generation of extra operators in 
(^-dimensions which do not have counterterms available from the original interactions. However, 
they differ from the generation of non-renormalizable terms in a theory by the fact that they 
exist only in (i-dimensions and vanish in the limit to two dimensions. Hence, they are accorded 
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the name of evanescent operators. As such they affect the extraction of the true renormalization 
group functions which is a problem we address in the three loop calculation. This absence of 
multiplicative renormalizability in Thirring like models has already been recognized in the work 
of jfq, 2C, 21, ^] upon which a lot of our calculation is founded. Though in the NATM was 



only discussed at the level of the two loop (3- function and anomalous dimension, whilst [21 



extended the wave function renormalization to three loops. These papers were important for our 
determination of the /^-function, anomalous dimension and mass anomalous dimension at three 
loops, which is our central aim. Due to these technical reasons, which will be explicitly detailed 
later, we concentrate on presenting a comprehensive discussion on the full renormalization of 
the NATM. Moreover, due to the nature of the renormalization we have also found it necessary 
to comprehensively study a variety of mass independent renormalization schemes, aside from 
the widely used MS scheme, to ensure consistency from other considerations, |l5| , ^l| , |22|] . This 
is guided very much by the fact that in two dimensions Thirring models are equivalent by, for 
instance Fierz rearrangements, to other four-fermi models such as the Gross Neveu model for 
certain values of Nf and N c . These properties need to be preserved in the quantum theory. 

The paper is organised as follows. The basic properties and relevant formalism is reviewed 
in section 2 in the context of the known renormalization of several two dimensional four-fermi 
models. The fundamental calculational details are discussed in section 3 prior to its application 
first to the Gross Neveu model with an SU(Nf) symmetry in section 4. Here we reproduce the 
known three loop MS renormalization group functions prior to presenting the detailed renor- 
malizion of the the NATM in the MS scheme in sections 5 and 6. We discuss the relation of 
these results to other renormalization schemes in section 7 before providing our conclusions in 
section 8. Various results relevant to the main discussion are provided in an appendix. 



2 Preliminaries. 

The form of the Lagrangian for the NATM which we will use is 

L natm = fyU^U _ m ^H^iI + | ^/T ; ^ iJ ) 2 (2.1) 

where we take Dirac fermions ifj l1 with an SU(Nf) internal symmetry, 1 < i < Nt, living in the 
fundamental representation of the colour group SU(N C ), 1 < I < N c . The generators of SU(N C ) 
are T a with 1 < a < (N^ — 1) satisfying the Lie algebra 

[T a ,T b ] = if abc T c (2.2) 

where f abc are the usual structure constants. With these definitions the fermion field is endowed 
with the same structure as quark fields in QCD. The more astute reader will have observed 
that the sign of the coupling constant, g, is opposite to the normal convention. One reason for 
this resides in the relation the interaction has with the Gross Neveu model which we discuss 



presently. Thus in its present form (|2.l| ) will not be asymptotically free. For the major part of 
the calculation the issue is one of performing the renormalization of the model and we will be 
able to restore the correct convention by the change g = — A at the end. Although the main 
focus will indeed be on the NATM, the calculation will run in parallel with a similar calculation 
in the SU(Nf) Gross Neveu model, [[!]]. There are various reasons for this. First, each possesses 
a four-fermi interaction with the difference residing in the presence of additional Lorentz and 
colour group structure in the former case. More concretely, its Lagrangian with the same fermion 
definitions is 

LS n = iftffl? ~ rnff + |(^V) 2 (2.3) 
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where, with this choice of coupling constant sign convention, the model is asymptotically free. 
Therefore, the integration rules for the Feynman integrals are the same in each case which will 
provide us with a useful cross-check on calculations in the NATM. In other words, since the 
d-dimensional values of the Gross Neveu graphs are known to three loops, [p!8| , ensuring that 
these are first reproduced exactly gives us confidence that the integration routines we use are 
correct. Second, for various values of the parameters of each model the theories are equivalent. 
These and other equivalences for the most general two dimensional four-fermi interaction have 



been summarized in [15|. Here we recall two which are relevant for this article. First, the abelian 



limit of ( [2.1D gives the abelian Thirring model, (ABTM), which has the Lagrangian 

L abtm = ^i^i _ + |(^*yy)2 . (2.4) 

When Nf = 1 one can show by Fierz rearrangement that 

(^V) 2 = - 2{^f . (2.5) 
Therefore, the Nf = 1 abelian Thirring model is classically equivalent to the Nf = 1 Gross 



Neveu model. Assuming that this is also valid in the quantum theory, |[^, 2C], means that the 
renormalization group functions have to be the same for Nf = 1. Moreover, ensuring that this 
is the case is tantamount to preserving the Ward identities which lurk within the quantum field 
theory. A second equivalence exists between the Nf = 3 Gross Neveu model and the Nf = 4 
non-abelian Thirring model which was discussed in fl5|| . Again we will assume that this is a 
property of the renormalization group functions. 

Indeed in this context it is worth recording the renormalization group functions as they stand 
prior to this work for each model. Although we have defined our model as having an SU(N C ) 
colour group, the 2-loop /3-function and anomalous dimensions are available for the model with 
a classical Lie group, 0, Q. In our notation, these are 

m = + nm^i + 0{gi) (26) 

and 

7(s) = _ wmmi + 0(g3) m 

where we recall that to this order the /3-function is scheme independent and the anomalous 
dimension, 7(5), is in the MS scheme. For reasons which we explain later, we will concentrate 
on the group SU (N c ) whose Casimirs are 

<h{R) = , C 2 (G) = N c , T(R) = \ . (2.8) 

However, the abelian limit is defined when the renormalization group functions are expressed in 
terms of the Casimirs of a general classical Lie group, by setting 

C 2 (R) -> 1 , C 2 (G) -> , T(R) -> 1 . (2.9) 

(For more insight into the relation between the abelian limit of a non-abelian theory, see |23|| ) . 
Hence, the abelian Thirring model renormalization group functions are quite simply 

P(9) = 0(g 4 ) , 7(5) = - M- + 0(g 3 ) . (2.10) 
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It is worth noting that there exist several all orders results for these renormalization group 



functions, [24]. These are 



A similar formula for the /3-function of the NATM when Nf = 1 has been deduced in |25f| . For 
the Gross Neveu model the renormalization group functions have been computed to two loops 
in g, [fi| and to three loops in MS in 0, H H] as 

m = _ + m + em*! + (2, 2 , 
7(9) = _ + w - m - ad + 0(/) ,, 13) 

Z7T 87T Z 327T' 5 



for an SU(Nf) internal symmetry. (We note that the results of [17, [Tg] corresponded to a Gross 
Neveu model with Majorana fermions in an 0{Nf) symmetry group. For an even value of this 
Nf one can complexify the fermions to ensure they lie in SU(Nf).) Setting N c = 4 and Nf = 1 



in (|2.6D and (2/7), we find 

x.'/) = — + 4 + °(s 4 ) , 7(5) = - ^4 - ()[ <-i ,) (2j5> 



and iVjp = 3 in (2.13), we have to two loops 

0(9) = - — + 4 + °(/) > 7(5) = - |4 + °(5 3 ) (2-16) 

which makes explicit the earlier equivalence after accounting for the relative minus sign in the 
coupling constant definitions. The apparent non-equality between both expressions for j(g) is 
accounted for by the realisation that we are comparing the anomalous dimensions in models 
with four and three fields respectively. 

It is worth commenting on the role of the mass term in each of ( |2.1| ) and ( |2.3[ ) . Strictly its 
presence breaks the continuous chiral symmetry possessed by each model, [||. Though in the 
case of the Gross Neveu model, which is the same as a model with an 0(2Nf) symmetry, it has a 
discrete chiral invariance. Therefore, the equivalences we discussed above will only be assumed 
to hold in the massless case. However, if one is considering perturbative calculations in this 
instance one inevitably will encounter infrared divergences. These can be regularized by using 



a modified or effective fermion propagator as discussed in [15] which is 

if 



(p 2 — m 



21 



(2.17) 



In this case m is not a parameter of the original theory and does not enter the corresponding 
renormalization group equation. Indeed the origin of m in the denominator is not from the 
Lagrangian but put by hand into the Feynman integrals which result from a purely massless 
propagator. For technical reasons concerned with checking the construction of the renormal- 
ization group functions in a consistent way for various schemes in this case, we will need the 
full renormalization group equation for the massive model. Therefore, we have performed the 
renormalization of each model in two separate cases. One which we refer to as the massless 



NATM, where the propagator ( 2.17 ) is used and another where the propagator is 

i(jj> + m) 
(p 2 — m 



2 . (2.18) 



5 



which we will refer to as the massive NATM. In this latter case m is treated as a parameter of 
the theory since fliPf ) follows from (|2T|) and it will appear in the full renormalization group 
equation. The Feynman diagrams will of course be the same for each case. 

The ultraviolet infinities will be regularized by using dimensional regularization. In this 
approach the dimension of spacetime is taken to be d = 2 — e where e is small. Integrals are 
performed in d-dimensions with the poles in e absorbed into the appropriate renormalization 
constants. However, in this analytic continuation the 7-algebra requires special treatment, 



pq , 27, 15, 21, |2^] which we briefly recall. In d-dimensions the set of products of 7-matrices 
defined by 

pMip-./xn = 7 [w 7 M2 _ _ _yxn] (2.19) 

where [. . .] denotes total antisymmetry, form a complete basis in the spinor space, |2(|. Clearly 
this space is infinite dimensional, only becoming finite when d is restricted to a positive integer. 
Therefore, in dimensional regularization of theories with fermions one, in principle, ought to 
decompose all 7-strings into this basis. For theories, such as QCD, the nature of the interaction 
and its ultraviolet structure ensures that the higher T( n ) do not arise, [21]. The situation 



for two dimensional four-fermi theories is different in that interactions of the form (ijjT/ n ^) 2 
will be generated in the renormalization. One consequence of this is that the theory is not 
multiplicatively renormalizable in principle and to extract the correct renormalization group 
functions requires care. We remark that this is a feature of dimensional regularization and if 
one were, for example, to use a cutoff in strictly two dimensions, the generation of these extra 
interactions would not arise. These are the evanescent operators we referred to earlier. In the 
context of the ABTM their presence has been noted before in (l5|, 21, ^] and the most general 



(abelian) scenario discussed in [15, ^0|. In both instances the analysis has been to two loops, 
though in [22| the consequences for the three loop wave function renormalization were also 
discussed. As an aside we remark that whilst the three loop Gross Neveu model renormalization 
has been performed, primarily to deduce the renormalization group functions at this order, the 



issue of whether evanescent operators are generated in this model had not been resolved, [18, GJJ. 
That they do not at two loops has already been established. Though it was noted in p^] that 
their generation at three loops only affects the form of the four loop renormalization group 
functions. As our goal is the full three loop renormalization of ([O]) we need to take account of 
the subtlety of the evanescent operators presence. Therefore, we recall the important features 
of the earlier explicit analysis of |]l^] . At the outset one can acknowledge the existence of such 
operators by including them in the most general possible Lagrangian, each with its own coupling. 
To make this more concrete we will focus for the moment on the ABTM before remarking on 
the NATM later. In this instance, the d-dimensional Lagrangian Q2.4|) now becomes 



1 00 2 
L abtm = ^i^i _ + -£ 5a (^rgpV) (2-20) 

" =0 



where all the fields and parameters are bare and g a are the couplings of all the possible four- 
fermi operators. For notational purposes in general discussion, we will use Roman letters to 
refer to only evanescent operators and distinguish the coupling of the original operator by g. 
Thus for the ABTM we take g% = g and setting g\~ = one recovers the original Lagrangian, 



(2.4). Moreover, when we consider the renormalization group functions for the model with an 
infinite number of couplings, we introduce the convention that if the argument of a function is 
the original coupling g, then it means that it corresponds to a renormalization group function 
evaluated in the case that g\~ = 0, for all k. Greek indices will correspond to both original and 
evanescent couplings and the letters n and k will be reserved for dummy indices in sums for the 
respective cases. 
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If, for the moment, we assume that the practical task of computing the Green's functions to 



some order has been achieved, then the renormalized parameters of (2.20) can be deduced with 
respect to some renormahzation scheme. So, for instance, the renormalized couplings will be 
a function of e and a finite number of other couplings. Therefore, the resulting /3-function for 
each coupling will depend on the others. Ordinarily in a theory with more than one coupling 
this is the standard method. However, since the couplings correspond to operators which do 
not exist in the limit d — > 2, one might naively believe that one can recover the renormahzation 
group functions of the theory by setting = 0. The resulting /3-function of the coupling g, 
however, will not correspond to the true /3-function of the theory. To understand this, it is 
best to consider the underlying renormahzation group equation of the general model, fl2.20| ). 
For a renormalized n-point Green's function G^ (p,m, /i, g K ) where p represents the external 
momenta and /U is a mass scale introduced to ensure the coupling constant remains dimensionless 
in d-dimensions, it satisfies, in our conventions, 



doc o f~. 

^ + 2_^W(^)^ + lm{g a )m— + -j(g 



dm 



G (n) (p,m,/i, 5/3 ) 







(2.21) 



where P( K )(g a ) are the naive /3-functions of the theory and 7(<? Q ) is the naive anomalous dimen- 
sion. In terms of the respective renormahzation constants, and Z/ K ^, these are defined to 
be 

d d 
P( K )(9a) = g K H—}nZ( K) , j(g a ) = fj,— ]nZ^ . (2.22) 



(Our convention for the 7(<?k) term of ( 2.21 ) is consistent with the definition of the perturbative 
functions, ( 2. 10| ) .) By considering the nature of Feynman diagrams, it is simple to observe that 
in the limit gk — ► 0, that the naive evanescent /3-functions do not vanish. Therefore there will be 
a remnant contribution in fl2.21|) in this limit, which would contribute to the true renormahzation 
group functions, which satisfy 



Q d d ti 



G^(p,m,^g) 







(2.23) 



where (p,m, /i, g) now corresponds to a finite n-point Green's function in strictly two di- 
mensions. In other words, the effect and contributions from the evanescent operators in higher 
loop corrections are not properly decoupled in the nai've /3-function, (3{g). Before reviewing how 
to do this we note our calculational strategy in light of the above remarks. Clearly there are 
two avenues of computation. The first of these is to take the most general Lagrangian with 
evanescent couplings truncated at the appropriate order which is dictated by the number of 
loops the renormahzation will be performed for. One then carries this out, computing the naive 
renormahzation group functions before proceeding to the extraction of the true renormahzation 
group functions. However, it is clear that this will necessitate a huge amount of calculation since 
gt 7^ 0. Since at the end of the computation the evanescent couplings will be unimportant, this 
provides us with the second approach which is the one we follow. We set g^ = at the outset. 
This reduces the amount of calculation from the point of view of Feynman diagrams but we 
will still be able to deduce the naive functions (3 K {g) and 7(5) which will now only depend on 
g. The procedure is to begin with the bare operator with a bare coupling g and perform the 
one loop renormahzation of the 2 and 4-point Green's functions. This will generate one or more 
evanescent operators whose coupling will be 0(g 2 ). To perform the two loop renormahzation, 
the Lagrangian with these new operators is used with the evanescent operators being treated 
as if they were vertices in the original theory. These newly generated vertices at this order are 
included in the subsequent Green's functions and so on. 
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The correct renormalization group functions in the chosen renormalization scheme are re- 



covered by applying the projection technique of 28] and discussed in the four-fermi context 



in [15, p0| ]. It is based on the observation that the insertion of an evanescent operator in a 
Green's function is not independent of the insertion of the relevant or physical operators in the 
same Green's function. More concretely, if we denote an operator by O k and its normal ordered 
version by J\f[0 K ], [p9fl , then within the context of a Green's function 

d d xM[O k ] = J d d x( K p^ k \g)M[ii)^ - m^V + 2gO{\ 

- pgXg)N[mW] + C^{g)N[0 1 ]) (2.24) 

/ gi=0 d=2 

where we take O k = ^(^T^ip) 2 and p^(g), Pm\g) and C^ k \g) are the general projection 
functions which will be computed order by order in perturbation theory. To do this one inserts an 
evanescent operator in either a 2 or 4-point Green's function and calculates it to the appropriate 
loop order. In general this Green's function with the operator insertion will be divergent, but 
one absorbs the regularized infinity into the composite operator renormalization constant which 
is available, with respect to the renormalization scheme being used. After renormalization one 
expresses this Green's function in two dimensions with = 0. The result will be non-zero. 
This exercise is repeated for the insertion of each of the operators of the right side of ( |2.24| ) and 
the finite value after renormalization is then used to determine the projection functions p^ k \g), 
Pm\g) and C^ k \g) perturbatively, [15]. The practicalities of this procedure will be clarified 



later when the explicit calculations are performed. Therefore, these projection functions are a 
measure of the presence of the evanescent operators in the appropriate Green's functions. If one 
now considers the action, S, of the theory (p.4|), 



S = [ d 2 xL ahtm (2.25) 



then it is simple to deduce, pq] , 

^dp~, = J d d x^j(g a )Af ,,;]),: - + 2 ^ g K K 

oo 

- 7m(Sa)A/-[mV^] + J2^)(9a)M[0 K ]\ . (2.26) 

K=0 / 

However, substituting for the relation of the evanescent operators to the physical ones, results 



in the simple relations, [15] 



M = /%) + T,C^(g)P {k) (g) 

oo 

7(9) = 7(5) + Y,P {k) (9)P { k)(g) 



k=0 

00 



1m (g) = lm{g) + (2-27) 
k=0 

from the realisation that the action of relevant operators obeys 

p,^L = J d d x( yl {g)M{ii>^ - + 2g0 1 ] - ^ m {g)M[m^] + P{g)M[0^j . (2.28) 

Following this procedure allows us to properly account for the evanescent operator problem. A 
similar reasoning, which we also found instructive in understanding the renormalization, arises 
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in the renormalization of the two dimensional Wess-Zumino-Witten model, 31]. There 
the evanescent operator which is generated occurs in the 2-point function and its presence can 
be properly accounted for in the construction of the renormalization group functions to three 
loops. The issue of evanescent operators in the context of renormalization has also been studied 
in 10211331. 



Although we have concentrated on the massive models in the above, one can easily recover the 
massless case by setting m = in the appropriate equations. Further, we have only considered 
the ABTM, since it is simpler than the non-abelian case. For the NATM, the procedure is 
exactly the same. The only difference being that we have to account for the presence of the 
generator T a . For a general colour group there appears to be no simple way of proceeding. 
Indeed a systematic study of group theory for general classical and exceptional Lie groups in the 



context of multiloop calculations has been provided in [ 34 ] which illustrates the deep complexity 
of the problem. Instead we follow the Cvitanovic procedure, [35|, and restrict ourselves to the 
colour group SU(N C ). For this group the generators satisfy 



IJ 1 KL 



1 



(2.29) 



which means that for each T( n ) there are two possible evanescent operators. Moreover, it will 
lead to a substantial simplification in the computation where we use ( |2.29 ) to decompose the 
interaction of ( |2.1[ ). After completing the calculation at each loop order, we then reconstruct 
the generator products. However, the full colour basis requires the inclusion of the identity I®I 
as well as T a T a for completeness. Therefore, we take as the generalized Lagrangian 



L 



1 00 r 2 

2 E 9,0 (^T^jj^) + g Kl [i, u T {K) Tfj^ J 



K = 



(2.30) 



where the second subscript in the evanescent coupling counts the number of generators. For the 
NATM, we now reserve g = gu for the coupling of the physical operator, with the projection 
functions now denoted by p^ kt \g), Pm\g) and C^ kt \g) with i = or 1. Although we have 
now explicitly specified SU (N c ) one can reconstruct the colour group Casimirs for an arbitrary 
classical Lie group at the end of the renormalization group function calculation by applying the 
results ( |2.8| ) where T(R) always accompanies one power of Nf. One issue which might arise with 
this is the appearance of Casimirs other than products of (|2.8| ) . From the nature of the Feynman 
diagrams, it is easy to convince oneself that new Casimirs, such as those which appeared at four 
loops in the MS /3-function and quark mass anomalous dimensions in QCD, || ||], will not occur 
at three loops in the NATM. 



3 Calculational details. 

We now turn to a discussion of the technical aspects of the algorithm used to carry out the full 
three loop renormalization of (|2.l| ) . We performed the calculation by a semi-automatic symbolic 
manipulation method using several languages and packages. For instance, the Feynman diagrams 
were generated using the package Qgraf, p6|], and then encoded in the language Form, p^j, 
as it was the most favourable package to perform the tedious amounts of algebra associated 
with the 7-matrices which we now discuss in detail. The advantage of using such packages 
can be fully understood from the nature of the problem we have already discussed. We are 
primarily interested in the structure of the Gross Neveu and NATM and to a lesser extent the 
ABTM. Our programmes are organised such that the graphs are encoded with a general vertex 
function carrying all the necessary Lorentz, spinor and group indices. One then specifies the 
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model by calling the routine which substitutes for the vertex peculiar to that model. This degree 
of flexibility has been fundamental in checking the results for the Gross Neveu model which is 
taken as the basic reference point. For any topology the evanescent vertex operator can also be 
implemented by using the same encoding for that topology but with the basic vertex replaced 
by the evanescent operator Feynman rule. This ensures that no graphs are accidentally omitted. 
For completeness we record the structure of the Feynman rules we have used for the vertices in 
each model in addition to (2.17, 2~T§| ). For the Gross Neveu model the basic interaction of (|2~ 
leads to the Feynman rule 



HJ 



^S^dJSJ - 5 il 5 kj 5j5j 3 



a "7 



(3.1) 



for external fields tp^, ify, i^ k and tjr s labelled anti-clockwise. Whilst the NATM has the Feynman 
rule for the set of fields, if) l J , $J ' , ip. 



,1 jJJ j,kK and 



>fj 



P > ^7 



I* Pry S 



and for the evanescent vertices 

ig [^5 kl 5 JI 5 LK T 

and 



(n)a 7 



6 ll 5 k ^ LI 5j K r {n) Jr^ 



The abelian limit of the original vertex is 



a i/i 7 



(3.2) 

(3.3) 
(3.4) 

(3.5) 



To appreciate the initial problems with the 7-matrices, we recall that each propagator has 
one 7-matrix whilst for the NATM each basic vertex contributes two more 7-matrices. Thus 
for any three loop topology there are at most fourteen 7-matrices. For the massive model each 
topology will also involve sums of less numbers of 7-matrices as well. From the fact that there are 
four free spinor indices the overall result will be the tensor product of two 7-strings. Given the 
nature of the vertex the two extremes of this tensor product will take the form • • • 7 
and 7 W 



M13 



> 7 



•7 



in 



■ j ur where the nature of the explicit contractions are not important for 
the present illustration. As we need to re-express the final value for the integral after integration 
in terms of the T( n )-basis we need to have a systematic way of reducing these strings to products 



of the form 7 Ml • • • 7' 



7/il ' ' ' lllr 



< 7. In the two extreme cases there are no free Lorentz 



indices on the strings. For the former it is clear that the string of thirteen 7-matrices will 
have six pairs of contracted indices. In the latter if there is a pair of contracted indices in one 
string then there will be a pair in the other string also contracted. Otherwise the contractions 
are across the tensor product. In general, rules for the treatment of the re n \-basis have been 



developed, [p7f| , and then applied to the three loop 2-point function, [21, g^[. We have not used 



the majority of these as they prove too cumbersome to implement symbolically. Instead we used 
a minimal number of their properties in addition to the (i-dimensional Clifford algebra 



{7/^,7^} = l<rf v . 

Our algorithm was based on ( |3.6| ) to shuffle contracted 7's together via, for example, 

7^7^ . . . 7^7^ = — 7 M 7 I/1 . . . 7^7^™ . . . + 2~f Unr y Ul . . . Y"' 1 
Then contracted pairs were removed with the usual simple identities 

7V7m = ~(d-2)Y 



(d-4) 7 v + ±rT 

- (d - 6)7V7 P - krf u l P + W - W 



(3.6) 

(3.7) 

(3.8) 
(3.9) 
(3.10) 
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This process was repeated until one was left with tensorial 7-strings which only had Lorentz 
contractions across the tensor product. At this point we decomposed each 7-string into its ]TV n ) 
basis, which requires inverting the simple definition, ( p. 19 ), This is achieved as follows by using 
the lemmas, [fill ^71, 



(n) i (n+1) ' i ' ' ' (n— 1) \ ' ' 

r=l 
n 

' (n) (n+1) / < v / ' (n— 1) ' .->._; 

r=l 

Beginning with the product 



rf 1)7 * (3.i3) 

(1 



we evaluate it in two ways. One way, of course, is simple in that = 7^ and thus the 
expression is a string of two ordinary 7-matrices, 7^7^. However, using the lemma once this is 
equivalent to 

r (S + ^ r (0) (3-14) 

where FVq) carries only spinor indices. It is easy to observe that this method is iterative. Be- 
ginning now with r^7^7°" and multiplying ( |3.14 ) by 7 " and applying ( 3.1 2| ) gives the relation 



between rj^ CT and r y^^ v ~f a and so on. Moreover, the iteration is trivial to implement symbolically 
and the output converted into the appropriate substitution routine. Once in this r( n )-basis any 
spinor trace which is left can be carried out by noting that 

t r r (n) = for n + (3.15) 

and trr(°) = 2. In practice, though, the performance is more efficient if we take traces at an 
earlier point. 

Although this completes the algorithm for performing the 7-algebra of the basic three loop 
graphs, the generation of the evanescent operators means that when they are substituted into 
Feynman diagrams the presence of r7 n ) for n 7^ 1 needs to be handled. First, one could decompose 
each T( n ) via ( |2.19| ), and then use the above algorithm to reconvert into the rV n ) basis after loop 
integration. However, this is not efficient, especially when one is dealing with the decomposition 
of T( 5 ) since it will involve 120 terms which with its tensor partner is beginning to generate a 
sizeable number of terms before one considers the number of Feynman diagrams for a topology 
or indeed the number of vertices in each. The more efficient way to proceed is to exploit another 
T( n ) identity which contains the first equation of (3.10) as a simple case. It is, JIB] , p7| , 



(n) (m) Mi-.-Mn J \ ' / (m) 



whereQ 

f(n,m) = (-\) nm (-l)< n - 1 )/ 2 \(l+u) d - m (l -u) m ] . (3.17) 

du n L J u=0 

At the outset any diagram involving the inclusion of an evanescent operator will be a mix of 
T( n ) and ordinary 7's. The latter, however, are decomposed into F( n )'s by the earlier algorithm 
before ( |3.16| ) is applied at the end. Having ensured that the 7-algebra can be handled in terms 
of the basis T( n ) we note that the most general set of 4-fermi operators covering the NATM is 

\ (fHuY^ 3 ) 2 , \(f^ { n)Thf J ) 2 ■ (3.18) 



*We note that this relation appears in both J27]] and plj . In there appears to be a minor typographical 
error in a phase factor in the first equations of (A. 1.2) and (A. 1.3). The correct factor, which we have used, ought 
to be Em and not Ek which can be readily checked by evaluating special cases of m and n. 
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One recovers this basis in practical terms by inverting ( 2.29 ) in conjunction with the spinor in- 



dices of the r7 n )'s. For instance, if one is dealing with the four point function (ipa^ J l3 ip!^ K ip lLS ), 



then 



$il5kj 



a (rc.) 7 



2TfjTZ K 



+ 



5ij5kl 



-p(n) /3-p 6 
a (n) 7 



(3.19) 



The other major part of the calculation is the integration of the Feynman diagrams which 



is based on the earlier papers of 17, 18], which we follow here. We will compute each of the 2- 
point and 4-point functions in both the Gross Neveu model and NATM (massive and massless). 
The former will be at non-zero external momentum whilst for the 4-point function we will set 
all external momenta to zero. This is another reason for the presence of a mass term which is 



necessary to have a scale in the integrals. In [18] the diagrams were computed in d-dimensions 
with the substitution d = 2 — e only being made at the end. We follow this source here primarily 
so that we can check that each of our integral programmes gives the correct result in the Gross 
Neveu model before applying to the NATM. There is one complication though. In the Gross 
Neveu model the vertex involves the spinor identity. Therefore within a topology the 7-matrices 
contracted with the momenta are adjacent. This simplified the Gross Neveu calculation as one 
did not need to consider tensor Feynman integrals. These, however, are necessary in the NATM 
and we have had to rework the results of |l8| to determine the tensorial nature of the integrals 
which arise for each topology. By way of illustration we recall that a typical integral is 



k^k u k a lPq x q 



Iklq {k 2 — m 2 ) 2 ((k — I) 2 — m 2 ){l 2 — m 2 )(q 2 — m 2 ) 2 
The most general form that can be consistent with the integral structure is 



(3.20) 



(3.21) 



However, by Lorentz symmetry it is simple to deduce that X = Y = Z. Hence the unknown 
amplitude is readily deduced by contraction with the resulting integral determined from the 



results in the appendix of [18]. It evaluates to 



I 2 + 



(d + 3) 2 



m 2 A(0) 



I 



Ad{d + 2) 



(3.22) 



Throughout the computation two basic integrals emerge which are present in ( |3.22 ). The first 
is the simple tadpole result 

I = i [ — 1 - (3.23) 



k (k 2 



which is divergent in two dimensions. The second term A(0) is the zero momentum value of the 
more general integral 

J(k 2 ) 



A(p) 



where 



J(p) 



k({k-p) 2 
1 



(3.24) 



(3.25) 



Ik {k 2 — m 2 ){{k — p) 2 

We leave A(0) unevaluated since it is connected to the issue of renormalizability. If we were 
performing the calculation of the 4-point function at non-zero external momenta the function 
A(p) as well as other finite momentum dependent functions would emerge. For the theory to 
be renormalizable such quantities must not appear in the full evaluation of the Green's function 
multiplied by poles in e. If they did then they would correspond to some non-local counterterm 
in the Lagrangian which would not fit the usual criteria for renormalizability. It occurs in a 
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subset of the two and three loop topologies. The appearance at two loops is as a finite term but 
it will contribute at three loops when multiplied by a one loop counterterm. In [|l7|, 18 1 their 
cancellation to give a finite contribution in the full expression for the 2 and 4-point functions 
was verified to ensure that the Gross Neveu model is renormalizable to three loops. We have 
checked this again within the context of our semi-automatic calculation before verifying that 
in the NATM and ABTM the sum of all A(0)/e pieces yield a finite contribution as it ought 
to, to preserve the usual renormalizability criteria at three loops. Since we will be considering 
non-minimal renormalization schemes in a later section we will require knowledge of the finite 
parts in addition to the poles in e. Therefore we have been careful in the construction of all the 
tensor integrals similar to ( [3.22 ) not to discard any (i-dependence. In other words even though 
say (d — 2) 2 I 2 is clearly finite in two dimensions we have not omitted it from any amplitude in 
the tensor decomposition. 






Figure 1: One and two loop topologies for the 2-point function. 



H & 

Figure 2: Three loop topologies for the 2-point function. 

Finally, we recall that the topologies for the various diagrams can be readily constructed. 
To three loops the basic ones for the 2-point function are given in figures 1 and 2. We have not 
labelled the graphs with arrows for the charge flow as to include all possibilites would increase 
the number of graphs significantly. Suffice to say that at one and two loops there are a total of 
24 graphs for the 2-point function and at three loops a total of 7 graphs. For the 4-point function 
the topologies are illustrated in figures 3, 4 and 5. In these cases there are more graphs due 
to the existence of various scattering channels which are not illustrated. In total there are 196 
one and two loop graphs and 213 three loop graphs. For both Green's function these figures for 
one and two loops take into account the contributions from diagrams with 2-point counterterms 
and evanescent vertices. We have not counted as separate those graphs generated by the vertex 
counterterm. 



4 Gross Neveu model. 



Before considering the renormalization of the NATM in detail we address the multiplicative 
renormalizability issue of the Gross Neveu model. In constructing our algorithm we have been 
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careful to check that the d-dimensional expression for each set of graphs contributing to the 
same topology agree with those given in [|l^]. Although that calculation was carried out for 
the model with 0(Nf) Majorana fermions, the observation that SU(Nf) = 0(2Nf) allows us to 
compare the results. Moreover, the comparison is made only for the <g> Tr^ type vertex. 
In [18 1 and [|1]] it was assumed that there were no other structures available which had not 
been justified. However, if new divergent structures did appear at three loops, they would not 
invalidate the three loop MS renormalization group functions computed in 17, [l^, [H^] • With the 
machinery of the r( n )-basis now available we have recomputed the full 4-point Green's function 
in the SU(Nj) Gross Neveu model. It turns out that at three loops several topologies give rise 
to non-zero contributions involving T( 3 ) g) T( 3 ) both for the massive and massless models. It 
could, of course, be the case that their sum is finite in two dimensions. However, we first record 
the values of the T( 3 ) £ED T( 3 ) part of the contributing graphs as 



G 4 i(0 3 ) 

G45(0 3 ) 



3 „4 



I*g 



, G 42 (C> 3 



3„4 



6d 



2d 



(3d 



G 44 (0 3 ) 
4)/V 



6d(d - 1) 



(3d - 4)lV 
12d(d- l) 2 

G 47 (0 3 ) 



3d(d - 1] 



(4.1) 



where the subscripts correspond to the graphs of figure 4. It is worth observing that these extra 







Figure 3: One and two loop topologies for the 4-point function. 

contributions in the T( 3 ) (g) r( 3 ) sector only arise in those diagrams with vertex subgraphs. It is 
easy to sum the contributions of the graphs to find the total is 



G (4) (Os) 



(3d-2)(d-2) 2 I 3 g 4 
Udjd - l) 2 



(4.2) 



which possesses a simple pole in e. Therefore, we need to include an extra set of evanescent 
operator counterterms in the original Lagrangian, (|2.3|) , which involved only bare fields and 
parameters. So, if the bare Lagrangian is 



- moVoV>o + \M%Y (4-3) 
then we introduce the Lagrangian involving renormalized fields and parameters as 

L§ n = iZ^W? ~ mZtZjptf + \g^Z g Zl(^^f + \g^Z^Z% (^T (3) ^) 2 (4.4) 
which will be valid to three loops and where we define the usual renormalization constants by 



■2 



m 



mZ n 



9o 



gZ g ^ . 



Tp = ipZ, 

We have checked that Z^, Z m and Z g are in agreement with the MS results of [17, 18 
Additionally, we have 



(4.5) 



>33 



48vr 3 f 



(4.6) 



14 



The form of Z33 is typical of the structure of renormalization constants associated with 
evanescent operators in that its leading term is not O(l). This can be appreciated by considering 
the Lagrangian analogous to ( p.20| ) where one begins with an infinite set of basis evanescent 
operators each with their own coupling. If one were to renormalize with that theory then the 
resulting renormalization constants would have a similar structure to Z g with the corrections 
involving all the couplings. However, in the limit — » 0, it is simple to see that the form 



would emerge. With (4.6) our observation is that the SU(Nf) Gross Neveu model is 
not multiplicatively renormalizable at three loops in dimensional regularization. Moreover, the 
lowest operators that can be generated in principle is Tr S ) ® Tr 3 \. This is due to the fact that if 
rVi) <S> or 1%) <S> Tr 2 ) appeared one would have to include counterterms for operators which 
do not vanish in the limit d — ► 2. This would be counter to the multiplicative renormalization 
in other strictly two dimensional regularizations. If we count the total possible number of 7- 
matrices which will occur in a 4-point function at successive loops up to three we find that these 
will be respectively 2, 4 and 6. It is clear therefore that the first appearance of a truly evanescent 
operator cannot occur before three loops which is what we find. 










Figure 4: Basic three loop topologies for the 4-point function. 



As will be discussed in more detail later for the NATM, the existence of such an operator 
only affects the renormalization of the model at the order after which it first appears. This 
is because the extra term in ( |4.4[) gives rise to a new vertex with a coupling involving Z33. 
Therefore the first diagrams where (|4.6j ) will be relevant are given in figure 6 where we denote 
the location of the 4-point evanescent operator by the cross in a circle. We note that one needs 
to take care when counting the g dependence of such graphs, since Z33 is 0(g 3 ). Whilst the 
anomalous dimension 7(g) has been computed to 4-loops in [38] in MS that result does not 
require the (g> rV 3 -j operator insertion in the 2-point function. It is easy to see this since the 
first Feynman integral, when computed, only contributes to the mass dimension. The integral 
involving the loop momentum in the numerator vanishes by Lorentz symmetry. Thus only the 
4-loop /3-function and mass anomalous dimension will require the computation of the graphs of 
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figure 6 in addition, of course, to the (large) set of diagrams arising from the original vertex. 
The detailed procedure to achieve this will be discussed in the context of the NATM where the 
same problem arises but at a lower loop order, [15]. 




Figure 5: Three loop topologies with tadpoles for the 4-point function. 




Z 33 

Figure 6: Contributions from evanescent operators to the four loop wave function, mass and 
coupling constant renormalization in the Gross Neveu model. 



5 Massive NATM in MS. 



As far as we aware there been two calculations dealing with the two loop renormalization of the 
(massless) NATM, [|14|, [L5} ]. The initial one performed in [14] used the auxiliary field formulation 
with a cutoff regularization which is not easy to extend to three loops. The latter one was 
detailed in the appendix of Eq] and used dimensional regularization. Using the elementary 
argument of the previous section to ascertain the form of the highest possible operator which 
could emerge from the 4-point function, it is elementary to observe that a T( 3 ) £ED T( 3 ) structure 
can emerge at one loop. This is evident in jy^] and as a preliminary to discussing the three 
loop calculation we will first reproduce the analysis of Jl5j in the context of this article. There 
will be several minor differences in our strategy. For instance, we will be interested in the mass 
renormalization which was not discussed before. Second, the results presented in [15] were for 
a general colour group and did not appear to appeal to the simplification manifested by ( |3.19| ) 
for SU(N C ). Although this is not necessary to two loops the problem of obtaining independent 
combinations of the colour generators is not trivial in the general case, (35]. Indeed it is tedious 
even at the two loop level. At the outset we note that any calculation we have performed at 
two loops has been checked with the NATM analysis of [15] including the coefficients of the 
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renormalization constants of the usual parameters and emerging evanescent operators, as well 
as the simple application of the projection formula to obtain the correct /3-function of 
Although we will extend the massless calculation of [15| to three loops, our strategy will be 
to first consider the massive model. The primary reason for doing this is to exploit non-trivial 
checks on the full renormalization group functions which become available in the massive model. 
Having understood the procedure for the massive model we proceed to the massless case where 
such a check does not exist. This is due to the different roles m plays in both models. In the 
latter case, since m is a regularizing parameter it does not appear in the full renormalization 
group equation. 

We begin by defining the renormalization constants we will compute, where the bare La- 
grangian is 



L 



natm 



Setting 



tpo = tpZ. 



m 



mZ„ 



go 



(5.1) 



(5.2) 



and allowing for all possible evanescent operators up to three loops, we will take the renormalized 
Lagrangian as 



L 



natm 



| 7 r 2 

k=0 L 



+ 



(5.3) 



where as in ( [4.4|) the Z k i are 0(g) except when k = 1 and i = 1. Although it would appear 
that there are thirteen terms up to level Trj\ ® ^(7) in addition to the original vertex, it turns 
out that only the operators involving rV 2r +i) f° r integer r, will occur, [pT| , p2j |, and the explicit 
calculation is consistent with this. To one loop the relevant diagrams are given in figures 1 and 
3. The tadpole, defined as a zero momentum insertion on a line, only contributes to the mass 
renormalization, and gives the value 



i(N 2 - l)dmlg 

For the 4-point function the sum of all channels gives 



(5.4) 



ig 2 i 



(A? - 1) 



N 2 

c 



d(d - 2)0, 



oo 



(N c 2 - 4) 



(N 2 

+ N c (d - 2)0 2 i + ' 



N 2 



<J 30 + 



d(d - 2)O i 
4) 



N c (3d - 2)0n 



o 31 



(5.5) 



Here, we adopt the practice that the Green's functions are written in terms of the (basis) 
operators O n o = ^(ipr^It/j) 2 and O n \ = ^(tl)T^T a tp) 2 whose definition does not contain any 
coupling constant dependence. By doing this one can easily define the counterterm appropriately. 



At this stage we choose to renormalize using the MS scheme and absorb only those parts of (5.4) 
and ( |5.5[ ) which are divergent. We have computed the coefficients of each operator as a function 
of d in order to consider other schemes which we will do later. Therefore due to the presence 
of the factors of (d — 2) in ( |5.5| ) counterterms are only required for On, O30 and O31 which 
have an odd number of 7-matrices. This is consistent with the earlier observation relating to 
chirality, [21, 22]. The generation of the latter operators affects the two loop renormalization. 
By counting the powers of the coupling constant in the Lagrangian, they need to be included 
with the two loop diagrams of figures 1 and 3, unlike in the Gross Neveu model where they 



17 



contribute at four loops. Therefore, in the 2-point function they are only relevant for the mass 
dimension through the tadpole of figure 6. For the 4-point function of figure 6, the insertion 
now represents the presence of both O30 and O31 . 

At two loops in the 2-point function there will be a contribution to the wave function renor- 
malization from the sunset diagram. Of all the graphs we consider, this, aside from its related 
graphs at three loops, represented the most tedious to evaluate. It was not possible to write the 
relevant integral in a closed form. Instead we treated it by expanding the integral in powers of 
p 2 around p 2 = and retained only those pieces which involve p* or mS a ^. These contributed 
respectively to the wave function and mass renormalization. In addition, we were careful to 
determine not only the pole parts but also the finite parts for the three loop calculation and in 
order to consider non-minimal renormalization schemes. We note that the wave function and 
other renormalization constants to three loops in the MS scheme have been collected together 
in appendix A. 

Computing the 4-point function now yields in addition to the operators O30 and O31, the 
new operators Oiq, O50 and O51. The latter two would be expected due to the simple 7-matrix 
counting, since the addition of a vertex introduces two 7-matrices from the new propagators 
and two 7-matrices from the extra vertex. For the one loop diagram of figure 6, the presence of 
T( 3 ) at the operator insertion also yields a maximum of ten 7-matrices which have the extreme 
decomposition of CST^ . The main issue with the two loop 4-point function is the observation 
that one cannot simply deduce the full /3-function from the naive renormalization constant Z\\. 
Ordinarily in the absence of the higher order operators the coefficients of the simple poles in e 
are simply related to the coefficients of the /3-function. In this case we would have 

N a 2 rfi 
P{9) = 0ii(g) = if- + (4iV c 3 iV) - 3iV c 4 + 12^-36)^-2 + 0( 5 4 ) . (5.6) 



As was pointed out in [15], this clearly contradicts the calculation of [14] since on general 
considerations the /3-function for a single coupling theory is scheme independent to two loops. 
The discrepancy between the two results is 3(i\, 4 — AN 2 + 12)g 3 / (16N 2 ) which is related to the 
/3-function of the operator T in the notation of where 

T = \^T a T h T c ^^T ahc -i^ (5.7) 

with 
and 

0A9) = ~ |J • (5-9) 



For SU{N C ) 



1 

SU(N C ) ~ 2 



(5.10) 



so that the second term is related to the spurious contribution to the naive /3-function. Clearly 
one needs a systematic method of accounting for the existence of such contributions and removing 
them from the true /3-function and the other renormalization group functions. 



This is provided for in the projection technique, |pq] , discussed in []T5| but which we recast in 
the basis of operators we have chosen. In addition to the naive /3-function, f3u(g), the evanescent 
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operators generated at one loop also have /3-functions, deduced in the usual fashion. To 0(g 3 ) 
these are 



030 fa) 



hi{g) 



(N 2 - l)g 2 

8N C 2 7T 

(K 2 - 4)g 2 

8N r 7T 



(5.11) 
(5.12) 



Although the general formalism was presented in section 2 its practical application requires 
explanation, especially as it will be applied to the three loop renormalization. As it stands, 
equation ( p. 24 ) represents the decomposition of the evanescent operator into a linear combination 
of the relevant original operators. The coefficients of this combination are perturbative functions 
of the physical coupling and the formula is only meaningful in the context of a Green's function. 
Substituting, say, the operator A/"[03i] into a 2-point function one evaluates the Green's function 
in two dimensions after any infinities have been removed. Then one inserts the right side of ( |2,24| ) 
into the same Green's function with the as yet undetermined projection functions and evaluates 
the Feynman diagrams to the same order before renormalizing and setting d = 2. However, it 
is easy to deduce that at leading order only the operator A/"[0n] will contribute for insertions 
in a 4-point function. Likewise Afliip^ip] and J\f[ipip] will be relevant for the leading order in the 
2-point function. This therefore provides the starting point of the perturbative iteration which 
determines p( kl \g), pm (g) and C^ kl \g), i = 0, 1. Indeed at leading order it is only the tree 
level insertions of the relevant operators on the right side of ( 2.24Q which will contribute. The 
insertion of these in the respective Green's function each evaluate to unity. 

Hence one needs only to compute the graphs of figure 6 for each operator O30 and O31 . For 
the 2-point function each is finite to this order and give 



^M[O 30 ]^) wf = 0(g) 
(i>N[O zl \i>) w f = 0(g) 



(5.13) 



and 



- + 0(g) 

IT 
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2M.7T 



+ 0(g) . 



(5.14) 



We have separated the contributions to the full 2-point function with the operator insertion into 
the piece involving j$ and the mass and denoted these respectively by the subscripts wf and 
to. However, for the 4-point function care must be taken to renormalize the operator insertion 
before setting d = 2. As this will become important later we record the insertion of each operator 
in the 4-point function to the 0(e) term 
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(5.16) 



where we have introduced the new mass scale fl 2 = 47re~ 7 /j 2 , with 7 denoting Euler's constant, 
to ensure that we are in the MS scheme as opposed to the MS scheme. The renormalization con- 
stants to remove the simple poles on the right side are provided by the fact that the expressions 
effectively represent the first stage in the renormalization of a composite operator in the normal 
fashion. There is a caveat here in that there is mixing into other basis operators. However, one 
can introduce a counterterm which is a linear combination of the basis elements. The upshot of 
this for higher order calculations is that when the renormalized operator is inserted one must 
also include the counterterms from the mixing operators. For instance, for O30 this will mean 
at the next order that a tadpole graph with insertions 4(2V? — l)Osi/(N c e) and — Osi/e have to 
be included. After this renormalization which we also take to be MS here, the contribution to 
the Green's function after setting the evanescent couplings to zero and d = 2 are, to this order, 



(^/W[o 30 ]#)| Oll 
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Hence, substituting these functions into the general formulas in (2.27) the actual renormalization 
group functions to two loops are 
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where the result of [14| is recovered. Moreover, we also have the MS mass anomalous dimension 
which has not been computed in this model before. 

We now turn to the details of the three loop calculations which builds on the above. We 
do not comment on the mundane task of the computation of the diagrams of figures 2, 4 and 
5 since their evaluation very much parallels any usual perturbative calculation in other models. 
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Instead we focus on the role of the evanescent operators. As the two loop calculation generated 
new operators their contributions must also be included. However, they will only enter through 
the topologies of figure 6 where the circle with a cross will now represent successively 0io, O50 
and O51. Additionally the counterterms generated for O30 and O31 are also included in those 
topologies as well as the original operator in the graphs of figures 7, 8 and 9. In figure 8 the 
insertion stands for the four possible combinations of O30 and O31. In these and the graphs 
with the original vertex we note that we have also included the usual wave function, mass and 




Figure 7: Graphs contributing to the NATM 2-point function from evanescent operators O^i in 
the two loop topology. 

vertex counterterms. Having detailed the extra evanescent operator contributions which have to 
be included it perhaps can be appreciated that such a calculation can only properly be attacked 
using symbolic and algebraic manipulation programmes. 




Figure 8: Graphs contributing to the NATM 4-point function from evanescent operators O^i in 
the one loop topology. 



It is not possible to include the values of each diagram in terms of the full operator basis. 
To illustrate the complexity of the result for one diagram, though, we record the structure of 
the three loop chain of figure 4. We found that for arbitrary d it is 
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Unlike the one loop result for the graph of figure 3 the evanescent operators which are even 
in the number of 7-matrices have, in principle, been generated with simple and double poles 
in e. However, it turns out that when all contributions are summed and evaluated that they 
have finite coefficients. Thus divergences are only associated with those operators which have 
an odd number of 7-matrices as expected from general arguments, pl| , [22| . Further, the new 
operators O70 and On appear which is consistent with our earlier argument. Renormalizing the 
contributions from the 2 and 4-point functions in MS determines the renormalization constants 












Figure 9: Graphs contributing to the NATM 4-point function from evanescent operators O^i in 
the two loop toplogies. 



given in appendix A. The resulting naive renormalization group functions for each operator can 
be deduced as 
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where aside from O30 and O^i we have only retained those terms of the evanescent operator 
/5-functions which are relevant for the true three loop renormalization group functions. To 
complete the calculation we need the subsequent terms in the projection functions. 

For the 2-point function this necessitates first of all computing the additional graphs to 
(ipAf[0 3 o]i>) and (tpj\f[0 31 }ip) as well as the new Green's functions (ipJ\f[O 10 ]'ip) , (i/jJ\T[O 50 ]iIj) 
and {ipN[0^i\ip) . As the details of the calculation of the latter three are very much akin to 
those of the former two at the previous order we merely record the results 
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For O^i their composite operator insertion into the 2-point function requires the computation of 
the graphs of figure 7 as well as the contribution from the graphs defined by the counterterms 
to the simple poles in ( |5.15| ) and ( 5.16| ) and their projection into the relevant operator. The 
upshot of including the additional graphs is that after MS renormalization we now have 
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To determine the extra terms of the projection formulas the insertion of the right side of ( 2.24 ) 
in a 2-point function must be included. Therefore, there is a contribution from the first graph of 
figure 10 after renormalization, from the A/*[On] term of ( |2.24j ). This enters with the coefficients 
C( 30 \g) and C^ 31 \g). Additionally the second graph, being 0(g), gives a potential contribution. 
However, it only contributes to the equation for pin(g)- Finally, the finite part of the right side 



N[i\|/3\|/1, N[\|/\|/] 





N[O u ] 



Figure 10: Graphs for the three loop projection formula in a 2-point function for the NATM. 
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Figure 11: Graphs for the three loop projection formula in a 4-point function for the NATM. 



The procedure for the 4-point function is similar. One also has extra graphs to consider 
involving the insertion of the relevant operators J\f[iip$ip] and Afltpip] and these are given in 
figure 11 for the right side of (2.24). Including counterterm contributions both for the evanescent 
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operator poles and the On vertex renormalization, the finite two dimensional Green's functions 
now become 
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Likewise, the sum of contributions to the right side give 
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for i = and 1. Hence, it is easy to deduce the projection functions to the order we require are 
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For the operators Oio, O50 and O51 we note that the one loop graphs of figure 6 with the 
appropriate operator insertions take simpler forms. Respectively, we have 
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The expression for (?/>?/\A/*[05i]V>V) is similar to that for O50 but involves more operators. As 
it does not contain any relevant new features we do not record it. Since the operator On does 
not appear on the right side of ( 5.43 ) and ( |5.44 ), then after renormalization and removal of the 
evanescent operators the contribution to the left side of the respective projection formula? is zero 
at this order. Therefore we have the simple results 

C W(g) = + 0(g 2 ) , C^(g) = + 0(g 2 ) , C^(g) = + 0(g 2 ) (5.45) 

which imply that at this order $io(g), P5o(g) and $5i(g) are not in fact needed for the renor- 
malization group functions. With the corrections to the other projection functions, it is now 
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possible to deduce the true MS renormalization group functions for the NATM at three loops. 
For SU(N C ) we find 
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It is worth noting that the term involving ln(/i /m ) has cancelled in the sum for 7 m (<?) 



Although we have detailed the projection technique of [15| for the massive NATM there is 
an alternative method of determining these renormalization group functions. In performing the 
calculation we were careful to compute the Feynman graphs as functions of d. Therefore, we can 
deduce the finite parts of the 2 and 4-point functions after renormalization. These must satisfy 
the full renormalization group equation ( |2.23| ) with the above functions. Therefore with 
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it is easy to verify that ( ^23|) holds. If we did not possess the MS functions then we could 
have solved for them by ensuring that the renormalization group equation ( 2.23j ) holds. For the 
massless model, since the mass is not a parameter of the theory it means that this avenue is 
not possible for us. This will leave us only the projection method which we have now shown to 
be consistent. In writing ( 5.50 ) we have retained only the contribution to the original relevant 
operator since after renormalization d is set to 2 and the evanescent operators removed. Finally, 
it is a simple task to rewrite the MS /3-function back in terms of the colour group Casimirs as 
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where we have also included the coupling A for comparison with the two loop results of | 



6 Massless NATM in MS. 

As the massive model clearly does not provide us with a fully chirally symmetric theory, we 
have also computed the MS renormalization constants for the model when the propagator fl2.17| ) 
is used. As the calculation runs very close to that of the previous section we concentrate on 
the new features. First, when one evaluates the contributing graphs to three loops with ( |2.17| ) 
the structure is simpler than the corresponding ones obtained by using fliDJ). Indeed the first 
graph of figure 4 now becomes 
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Comparing this with the expression for the same graph in the massive case we observe its decom- 
position into the operator basis involves only those operators with r^r+i) f° r integer r. This is a 
result of the chiral symmetry which underlies the choice ( 2.17 ). Moreover, this is also a property 
of all the other graphs to this order, so that the potential need for renormalization constants 
Z n i for n even and i = or 1 is automatically excluded. Completing the MS renormalization of 
the model as previously it turns out that the renormalization constants are the same as for the 
massive model with exception of Z m which, of course, does not arise at the outset. 

The computation of the full 'y(g) and (3(g) in the massless NATM now follows the use of the 
projection formula for the massive model whose practical use we have just demonstrated. Since 
m is no longer a true parameter of the model, the projection formula (2.24) clearly takes the 



simpler form, [15], 
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The absence of Pm (g) and, in principle, different values for the contributing terms to the 
Green's function structure might mean that the true MS renormalization group functions will 
differ from those of the massive model. However, as argued in Jl5|, 22] one would expect them 



to be equivalent. Therefore, it is important that this is checked explicitly. Though we need only 
do this for (3(g), as clearly j(g) will be unchanged. 

For the 4-point function there are again no contributions to the relevant operator in the 
Green's function involving separate insertions of 0io, O50 and O51 leaving us to concentrate on 
030 and 03i. For these, in order to illustrate the differences in the structure of the Green's 
functions which are needed to deduce C^(g) and C^(g), we note that their one loop forms 
are 
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which can be compared with those of the massive NATM calculation. Further, including all 
relevant counterterms the restriction of these Green's functions in d = 2 at two loops gives 
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which differ from those of the massive model. To deduce the final values of the projection 
functions, we record that now the right side of (|6.2| ) is simpler than ( p. 24 ) due to the absence 
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of p m (<?)• Explicitly, we have 
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for i = and 1. However, when solving perturbatively for the projection functions it turns out 
that their values to this order are the same as before. Therefore the MS /3-function is identical 
to that of the massive model as expected. The discrepancy between each set of formulas rests 
in the hidden role of the mass term in various Green's functions. 











Figure 12: Graphs contributing to the renormalization of [ipip] in the massless NATM. 



Although it may seem that for the massless model that there is no renormalization group 
function for the mass, it is possible to compute the anomalous dimension of the composite 
operator whose coupling would correspond to a mass. Clearly this is [tpip]. As this operator 
will be renormalized, one can define an anomalous dimension 7^(5) from its renormalization 
constant which can then be regarded as the mass dimension. Indeed the computation of 
the three loop quark mass dimension in MS in QCD was performed in the massless theory by 
renormalizing the 2-point quark Green's function with a [tptp] insertion, [13]. We have computed 
the relevant diagrams to three loops for both the massless Gross Neveu model, in order to check 
our integration routines, and the massless NATM. These are illustrated in figure 12 where we 
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have only shown those diagrams which are not trivially zero. Also, the cross indicates the location 
of the insertion of [ipip]. Therefore, if the insertion in the two loop tower of bubbles had not been 
on the top loop but the lower one, then the diagram is simply zero because in the massless model 
the one loop tadpole vanishes by Lorentz symmetry. Therefore if a topology has one tadpole then 
for a non-zero contribution its [ipip] insertion must be in that tadpole. Although it may appear 
that we have to introduce a new set of integrals to carry out the calculation, we have exploited 
the fact that since we are only interested in the [ipip] renormalization, the graphs need only be 
computed at zero external momentum. In this case it turns out that tensor integrals for each 
graph can be paired with a graph used in the evaluation of the 4-point function. The topology 
of the related graph is deduced by regarding the insertion as the location of two external legs. 
So, for example, the last graph of figure 12 is paired with the penultimate graph of figure 4. 
Similarly the first graph of the bottom row of figure 12 is related to the fifth graph of figure 4. 
This observation ensures that we do not have to unnecessarily undertake extra calculation. It 
turns out that performing the MS renormalization of these graphs that is identical to Z m of 
the massive model for both the Gross Neveu model and the NATM. To complete the evaluation 
of the associated true renormalization group function, we need to extend the projection formula 
of ( |6.2j ) to the case of bona fide operator insertions. In this case we take it to be 
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Here the last term has the insertion of unity since we regard this operator equation as being 
inserted in a Green's function of the structure (ip . [ipip]ip) which includes the original operator 
insertion itself represented by the dot. Therefore, one needs to be careful in distinguishing 







Figure 13: Graphs contributing to the massless NATM 2-point function with 
the evanescent operators O-^. 



b] insertion from 



which operator is which when considering the graphs contributing to the Green's function. To 
this end we have illustrated the relevant non-zero contributions to the two loop Green's function 
(ipAf[0}[ipip]ip) in figure 13. As we have detailed similar calculations earlier, we record the values 
of the Green's function for the operators Af[Oso] and A/"[C?3i] as those of other operators will 
only involve one loop graphs. We found, 



(ipN[o m ][ipip]ip) 
{il>Af[o 31 ][w$) 
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In 



m 2 



+ G(g 2 ) 
(6.8) 



where we note that the overall sign difference between these and the analogous results in the 
massive model is due to the fact that here we are inserting the operator [ipip] and not the negative 
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of this. Moreover, 
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+ 0(g) 
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(6.10) 



so that the numerical values of p^\g) and Pm°\g) differ. For the right side of the projection 
formula we have displayed the one loop corrections which need to be computed in figure 14. 
With their values the respective right sides of each projection formulas are 
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for O31, from which it is easy to deduce 
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Therefore, the MS 7^(5) is identical to 7 m (<?) of the massive model. 




N[iy3y] g 




N[O u ] 



Figure 14: Graphs for the three loop projection formula in a 2-point function with 
for the massless NATM. 



insertion 



7 Symmetric schemes. 

Although we have produced the renormalization group functions for both massive and massless 
versions of the NATM in MS there still remains the problem of reconciling the results with the 
equivalences with other models discussed earlier. The MS expressions do not represent results 
which are consistent with symmetries of the original model and the problem can effectively be 
traced to a hidden 7 s issue, [15, 22]. In other words, when one deals with models with 



continuous chiral symmetries in dimensional regularization this symmetry is not preserved in 
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MS. For example, in QCD one must treat 7 5 carefully in renormalizing quark bilinear currents 
involving 7 s . There the resolution is to use a renormalization scheme which does preserve 
the symmetry. This is achieved by first performing an MS calculation to the requisite order 
and then introduce an additional finite renormalization to ensure that in that case the Green's 
functions obey the (quantum) Ward identities, including the axial anomaly, |^9|. For the present 
calculation the general approach is the same. Though there will be several minor differences. 
One rests in the fact that the analogous finite renormalization is chosen so that one ensures 
the renormalization group functions in the first instance and not the Green's functions are the 
same in the equivalences between the various models. Further, in the present context, similar 
finite renormalization schemes have been established for the three loop anomalous dimension 
ensuring agreement of the equivalences between the ABTM and the Gross Neveu model for Nf 
= 1, |p2j . In essence that scheme was specified by removing in addition to the simple poles, the 
(numerical) finite part of the Green's function into the renormalization constant, p2]| . We stress 
that terms involving ln(/i 2 /m 2 ) were not removed so that one is still within the context of mass 



independent renormalization schemes. In [15, 22] this was referred to as the symmetric scheme. 



To appreciate the effect such a choice of scheme can have on the renormalization group functions 
we have determined them for the Gross Neveu model. To three loops for the two dimensional 
massive model, we found that 

m = _ w^M! + - m-w-w + ofe5) (7 ,) 



7(9) = - ' V, - - ' ~ "T " + °<f 4 ) < 7 - 2 > 



im(g) 



tt 2vr 2 8vr 3 

(2 N £ - l)g 2 _ (2N f - l)(2N f - 3)g 3 

8vr 2 8vr 3 
(2N f - l)g {2N f -l)(AN f -r)g 2 

2vr 8tt 2 
(8iY/ - 42N f + 37)(2iVf - l)g 3 



16vr 3 



+ 0(ff 4 ) • (7.3) 



Clearly, one now has a cubic polynomial in Nf at three loops in j m (g) in contrast to the quadratic 
that appears in the MS scheme. We can readily produce this result since we have been careful in 
the contruction of the integration routines to keep all the (/-dependence. The aim now is to extend 
our calculation to determine the scheme or schemes which preserve the known equivalences of 
the NATM. We follow the approach of 15, 22] and reconcile the results with the MS forms 



of the 0(2Nf) Gross Neveu model which possesses a discrete 7 s chiral symmetry. Though we 
stress that one could choose another scheme as reference. Further, since the choice of scheme is 
very much arbitrary we will take a general approach and introduce a parametrization. This will 
involve a set of variables which will be constrained by the equivalences. To illustrate this we will 
consider the massless model only and define the finite terms of the wave function renormalization 
to be 

4^M+^ + f). <M) 

Likewise for the coupling constant renormalization constant, ZnZ^, we take the finite part to 
be 

b n N f + b 12 N c + ^) ~ + {b 21 Nf + b 22 N f N c + ^ + b 24 N c 2 + b 25 + ^ . (7.5) 

An 0(g 3 ) term will affect the four loop form of the renormalization group functions. The form of 
the additional terms has been chosen to be consistent with the colour Casimir structure that is 
allowed. Therefore, the term of ( |7.4| ) is proportional to C 2 (R). Of course, these additional terms 
will affect the form of poles in e of the renormalization constants which we have not recorded. 
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In the context of the usual multiplicatively renormalizable theories the treatment of a scheme 

1.) Once one has 



change is straightforward and well documented. (See, for example, [40, 21 



one scheme, the renormalization group functions in another are determined by the relations 



7(5) = lid) 







20(g)— ]n£(9) 



(7.6) 



where 



(7.7) 



and the bar denotes the corresponding quantity in another scheme. In general the relation 
between the couplings in both schemes is determined by comparing the renormalized 4-point 
Green's functions using ( fT7^ ) after the function is deduced at the appropriate order in the 
2-point function. For the non-multiplicatively renormalizable NATM the situation is the same, 
though of course one has to first compute the new Green's functions in the scheme involving 
( |7.4| ) and ( |7.5| ). To understand fully the proper procedure for this we have first performed the 
renormalization in the massive model and then determined the renormalization group functions 
by finding the function £(g) before applying the formulae, ( |7.6| ). These could then be checked 
against the renormalization group functions deduced by ensuring that the renormalization group 
equation itself was satisfied. Repeating this first method for the massless model we have the 
general functions 
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which therefore provide us with the general scheme results 
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The MS results are of course recovered by setting the various constants to zero. 
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The two equivalences we wish to preserve are those discussed in section 2. For the one 
involving the ABTM and the Gross Neveu model at Nt = 1, we first need to convert the SU(N C ) 
dependence into the general colour group Casimirs and then apply (2.9). We find 



7(0) 



2vr 2 



N? + 4bnN? - 4bi 3 N f 



2vr 3 
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(7.12) 



where the latter two equations arise from choosing to remove the numerical part of the 2-point 
and 4-point functions respectively at three loops, similar to p2]]. With the caveat we had 
previously on the form of j(g) for the second equivalence, we find that the constraints in this 
case are 



= I6&12 + 5&i3 + 80w 2 i + 256u; 2 2 
= 81926f 2 + 51206i2&i3 + 8006 2 3 - 20486 22 
+ 480&26 + 4800^21 + 15360u>22 + 231 . 
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24 
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Clearly not all the general coefficients can be determined by these restrictions and the undeter- 
mined coefficients represent the remaining freedom in the choice of scheme that is possible. For 
instance, if we set 



15 

W21 = bu = 613 = , 622 = ^ ' &23 = 4 ' 
then the renormalization group functions will take the form 
(iV c 2 - l)N f f (iV c 2 - l)N f f 
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(7.15) 



In the abelian limit, we obtain (5(g) = and 7(3) = — Nf g 2 / (2ir 2 ) + 0(g 4 ) consistent with 
(2.14). Moreover, comparing j(g) with the choice made in p2| one observes that there are other 



schemes aside from the symmetric one of |22j for which the general form can be reconciled. 



Of course, choices other than (7.14) can be made. Finally, this completes our analysis of the 
basic renormalization group functions for the original massless non-abelian Thirring model. 



8 Discussion. 



We conclude by stressing that we have now completed the first three loop renormalization of the 
NATM in a dimensional regularization. Whilst it involved the careful extraction and treatment 
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of evanescent operators in the four point function we were able to determine the renormaliza- 
tion group functions of the model at a new order in perturbation theory in a variety of mass 
independent renormalization schemes. Moreover, from a calculational point of view we have 



also demonstrated how one interprets and applies the projection formulae of |15| , 29] beyond the 
leading perturbative order which was used to determine the 2-loop (5- functions of |lf|] . Although 
we have concentrated on the renormalization of the NATM, another important property we have 
established is the breakdown of multiplicative renormalizability in the Gross Neveu model at 
three loops in dimensional regularization. For the NATM this, of course, occurs at one loop, ijlSj] . 
In |2^] the consequences of the perturbative non-multiplicative renormalizability of the Gross 
Neveu model were discussed from the point of view of the (non-perturbative) large A^ expan- 
sion and the determination of critical exponents at a fixed point of the d-dimensional /3-function. 
Indeed it is worth reviewing some of these results now that the property has been established. 
Such exponents are fundamental, being scheme independent, and relate to experimentally mea- 
surable quantities. One ordinarily determines them from the (i-dimensional renormalization 
group functions through the e-expansion, using perhaps resummation techniques to improve the 
series convergence for applications in, say, three dimensions. The large expansion comple- 
ments such an approach in that perturbation theory is reordered so that bubble chain graphs are 
summed first. Indeed in the context of the Gross Neveu model the large Nf expansion has been 
successful in c?-dimensions in determining the wave function exponent r\ to 0(1/A^), [41, fj3|| , 



and the mass and /3-function exponents to 0(1/A^ 2 ), 45]. For the three loop perturbative 
MS results which had been computed there is exact agreement in the region of overlap since 
the exponents are related to the renormalization group functions through, for example, rj = 
7(3*) where g* is the value of the critical coupling in (^-dimensions. However, in light of the 
perturbative non-multiplicative renormalizability of the Gross Neveu model, it is worth first of 
all recording the basis of the large Nf formulation. Rather than using the Lagrangian (^1^) the 
starting point is the massless Lagrangian written in terms of an auxiliary field, a, 

2 

LS n = i&frb* + - — . (8.1) 

2ff 



One then develops the large Nf formalism using ( |8.l|) . With this form it does not appear 
possible to generate another operator of the form a iiVPr i\) l T] li {,p'^ % which would be analogous 



to the evanescent operator generated in Q4.4j ). For instance, if one considers the graphs of 
figure 4 which give rise to the evanescent operator, it is evident that a rV 3 ) (g) rV 3 ) structure 
cannot emerge in the 1/A^- expansion. With this apparent discrepancy between the renormalized 
perturbative theory and large Nf approach, it has been suggested in |^] that the simple relation 
V = lid*) needs to be modified in the perturbative approach to account for the presence of the 
evanescent operators in d-dimensions. In other words in perturbation theory one assumes that 
in (i-dimensions the Gross Neveu model possesses an infinite set of couplings g K for each operator 
\{"4>T ( K )ip) 2 corresponding to a Lagrangian which is multiplicatively renormalizable. Then the 
fixed point of the (i-dimensional theory is defined by /3( K )(<7*) = for all k. Whilst in practical 
terms this may be difficult to determine, if we assume that there is a solution then one can define 
the actual exponents by extending the naive relation for the wave function to r/ = 7(5^), ] p2j] , 
where now the renormalization group functions will depend on all the couplings. Turning away 
from this ideal scenario one can approach this infinite set of fixed points in an appropriate way, 



[22], by assuming that it is close to the one defined by g$ ^ and g* = 0. In other words the 
case where all the evanescent couplings are set to zero, |22]| . Several issues immediately arise. 
First, it is not clear whether the naive (3- functions of the (i-dimensional theory are the ones 
which define the fixed points. Second, the practical task of computing the corrections to the 
evanescent fixed points needs to be resolved. In reviewing these issues we believe the simplest 
strategy to tackle this problem is to perform explicit calculations. As the evanescent operators 
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in the Gross Neveu model do not become evident before three loops, the first stage would be to 
determine the four loop MS renormalization group functions. Second, since the renormalization 
constant for the evanescent operator is Nf independent its effect in the 1/Nf expansion will not 
become apparent at least before 0{X/Nf). In fact by computing the graphs of figure 6 it turns 
out that they are Nf independent. This follows from the fact that to obtain any Nf dependence 
one must have a closed fermion loop. However, from the nature of the new evanescent vertex 
such a loop in any of the graphs is associated with tr which vanishes. So for the /3-function 
and the mass dimension one can not access the contribution from this operator in large Nf prior 
to 0(1/Nf). Therefore to reconcile higher order large Nf results for critical exponents with 
explicit perturbative calculations will require a huge amount of calculation. 
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A MS renormalization constants. 



In this appendix we record the values for all the renormalization constants in the MS scheme to 
three loops in the non-abelian Thirring model. From the 2-point function, we find 
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For completeness, we also provide the values for the same renormalization constants in the 
abelian Thirring model in the MS-scheme. These were computed separately from those of the 
NATM, but agree with them when the abelian limit of ( |2.9| ) is taken. We have 
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From these renormalization constants, it is straightforward to compute the naive renormalization 
group functions to 0(g 3 ). They also agree with the abelian limit of the analogous functions in 
the NATM, apart from w (g). We find 
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Finally, we note that the values of the MS renormalization group functions to three loops are 
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